MIT-CTP 4103 



Is a color superconductor topological? 



Center for Theoretical Physics, 



Yusuke Nishida 
Massachusetts Institute of Technology, 
(Dated: January 2010) 



Cambridge, Massachusetts 02139, USA 



o 

< 



> 

o 
o 



X 



A fully gapped state of matter, whether insulator or superconductor, can be asked if it is topo- 
logically trivial or nontrivial. Here we investigate topological properties of superconducting Dirac 
fermions in 3D having a color superconductor as an application. In the chiral limit, when the pair- 
ing gap is parity even, the right-handed and left-handed sectors of the free space Hamiltonian have 
nontrivial topological charges with opposite signs. Accordingly, a vortex line in the superconduc- 
tor supports localized gapless right-handed and left-handed fermions with the dispersion relations 
E = lizvpz {v is a parameter dependent velocity) and thus propagating in opposite directions along 
the vortex line. However, the presence of the fermion mass immediately opens up a mass gap for 
such localized fermions and the dispersion relations become E — ±w\/m2~f"pf. When the pairing 
gap is parity odd, the situation is qualitatively different. The right-handed and left-handed sectors of 
the free space Hamiltonian in the chiral limit have nontrivial topological charges with the same sign 
and therefore the presence of the small fermion mass does not open up a mass gap for the fermions 
localized around the vortex line. When the fermion mass is increased further, there is a topological 
phase transition at m = \//?"+"A2 and the localized gapless fermions disappear. We also elucidate 
the existence of gapless surface fermions localized at a boundary when two phases with different 
topological charges are connected. A part of our results is relevant to the color superconductivity 
of quarks. 



PACS numbers: 21.65.Qr, 12.38.-t 



INTRODUCTION 



A fully gapped state of matter can be asked if it is 
topologically trivial or nontrivial. The most well-known 
topological state of matter is the quantum Hall effect in 
two dimensions (2D) [1-3], in which the time reversal 
symmetry is explicitly broken. Recently the time rever- 
sal symmetric extensions of the quantum Hall effect have 
been theoretically proposed and experimentally observed 
both in 2D [4-6][7] and 3D [8-10] [11-15], which are re- 
ferred to as topological insulators (the 2D topological in- 
sulator is also known as the quantum spin Hall effect). 
Generally the topological state of matter is characterized 
by the nontrivial topological charge of the single-particle 
Hamiltonian and accompanied by topologically protected 
gapless edge/surface states with linear dispersions (Dirac 
fermions). 

Another class of topological states of matter arises in 
superconductors. Although the ordinary nonrelativistic 
s-wave superconductor is topologically trivial, the weakly 
paired phase of the Px+iPy superconductor in 2D is topo- 
logically nontrivial [16] and Sr2Ru04 is its candidate ma- 
terial [17]. In addition to the gapless edge/surface states, 
the vortex in the topological superconductor supports 
gapless fermions localized around the vortex core. It is 
also known from the pioneering work by Jackiw and Rossi 
that the relativistic s-wave superconductor in 2D has the 
similar properties [18]. Such a system has recently re- 
ceived renewed interest because it can be realized on the 
surface of the 3D topological insulator in contact with 
the s-wave superconductor [19]. It is pointed out that 
the Balian-Werthamer state realized in the B phase of 
the superfluid ^He is also topological [20-23]. 

Such progress on the discoveries of the topological in- 



sulators and superconductors motivates us to ask the fol- 
lowing question: Is a color superconductor topological? 
In order to shed light on this question, we investigate in 
this paper the topological properties of superconducting 
Dirac fermions in 3D. In Sec. II, we start with the mean- 
field model Hamiltonian of the color superconductivity of 
quarks and point it out that the free space Hamiltonian 
can be characterized by a Z valued topological charge. 
In Sec. Ill, we compute the topological charge of the free 
space Hamiltonian when the pairing gap is parity even as 
a function of the fermion mass and show that its value is 
closely linked to the existence of gapless fermions local- 
ized around a vortex line. The low-energy spectrum of 
such fermions is also determined in this section. 



In Sec. IV, we turn to the case where the pairing gap 
is parity odd. By studying the topological charge of the 
free space Hamiltonian, we show that there is a topolog- 
ical phase transition as a function of the fermion mass, 
which is also reflected in the existence of gapless fermions 
localized around a vortex line. In Sec. V, we elucidate the 
existence of gapless surface fermions localized at a bound- 
ary when two phases with different topological charges 
are connected. Finally, Sec. VI is devoted to the sum- 
mary of this paper and implications of our results for the 
color superconductor are discussed. 
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II. PREPARATIONS 

A. Model Hamiltonian for the color 
superconductor 

We start with the following mean-field model Hamil- 
tonian for the color superconductivity of quarks [24]: 



Hcsc 



dx 



tplj {-ia ■ d + lim - pi) 6abSfg V'h,; 



(1) 

Here C = 17^7° is the charge conjugation matrix and 
wc assumed the same mass m and the same chemical 
potential n for all three colors (a, b) and all three flavors 
{f,g)- ^abjg is the pairing gap in the Lorentz singlet 
and even parity (J^ = 0+) channel and its color and 
flavor structure is specified by 



no supercurrent), one can choose A(a;) to be a real func- 
tion.^ Then T-L has the time reversal symmetry: 



with 







/ = I 7I73 



(5) 



Because H$ = leads to H{T^*) = E{r^*) and 
$'''T$* = 0, the spectrum is at least doubly degenerate. 
Therefore, we find that our Hamiltonian T-L belongs to the 
symmetry class Dili in the terminology of Refs. [25, 26].^ 
In Sec. IV, we will also considcir the case where the 
pairing takes place in the odd parity channel (J^ = 0~). 
Such a case is described by simply replacing A(a7) in 
Eqs. (1) and (3) with 7'''A(a;). The resulting Hamilto- 
nian has the same charge conjugation and time reversal 
symmetries and thus belongs to the symmetry class DIII 
again. 



Topological charge for class DIII Hamiltonians 
in 3D 



^abjg{x) = ^ Ai{x)eiabeifg. (2) 
i=l,2,3 

Ai = A2 = A3 7^ corresponds to the fully gapped 
color-fiavor-locked (CFL) phase and Ai = A2 = 0, A3 7^ 
corresponds to the two fiavor pairing phase where only 
four quarks are gapped. In either case, an appropriate 
transformation of ipaj by a real and orthogonal matrix 
in the color and fiavor space can bring the Hamiltonian 
(1) into the decoupled form i?csc = where 



ipj {—iot ■ d + 13m • 

1 



jdx 



\ J dx (V.t -^JC7^) 



-ia ■ d + (3m 

A*ix) 

\ Jdx-^tHj^j. 



ia - d — 13m + n 



^3 



(3) 



For example, in the CFL phase, we have Ai^2,3 = 

— A4 8 = ^Ag. Below we concentrate on one sector 

with the nonzero gap Aj ^ and suppress the index of 
j- 

The single-particle Hamiltonian H in the Nambu- 
Gor'kov representation has the charge conjugation sym- 
metry: 



According to Ref. [20], fully gapped 3D Hamiltonians 
belonging to the symmetry class DIII can be classified 
by a Z valued topological charge. The topological charge 
is defined for the free space Hamiltonian where A (a;) = 
A*{x) = Ao is a constant. Suppose the single-particle 
Hamiltonian in the momentum space Hp = e"'^''" H e'^ '" 
is diagonalized as 



Hp — Up 



Dp 






-D, 



(6) 



where Up is a unitary matrix and Dp is a diagonal ma- 
trix with positive elements. For the fully gapped Hamil- 
tonian, we can adiabatically deform Dp into the identity 
matrix, which continuously deforms the Hamiltonian into 
a "new Hamiltonian" Qp defined by 



(7) 



with properties Qp = Qj, and Qp = 1. Because of the 

charge conjugation and time reversal symmetries, a uni- 
tary transformation can bring the Hermitian matrix Qp 
into a block off-diagonal form [27]: 



Qp 



qp 
4 



(8) 



Qp = 1 leads to qpq^ = 1 and thus qp is a unitary matrix 
[U(4) in our case] . The topological charge is provided by 



c-^nc = -n* 



with 



C = 



-C7'^ 

C75 



(4) 



Because H^ = leads to n{C^*) = -E{C^*), the 
spectrum is symmetric under E o —E. Furthermore, 
when the phase of A(a;) is uniform over the space (i.e. 



^ If we did not make this choice, we need to modify the definition 

of r. 

^ Class DIII Hamiltonians are meant to have the charge conjuga- 
tion and time reversal symmetries with properties = C and 

= -r [20]. 
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the winding number of qp associated with the homotopy 
group TTz[U{n > 2)] = Z: 



N 



1 



247r2 



j dpe'i^i:x[{q-^diq^){q-'djqj,){q-^dkqp)] • 

(9) 

When a given Hamiltonian has a nonzero topological 
charge, such a system is said to be topological. Because 
Hamiltonians having different topological charges can not 
be continuously deformed into each other without clos- 
ing energy gaps in their spectrum, the topological charge 
defined in Eq. (9) classifies 3D Hamiltonians belonging 
to the symmetry class Dili [20]. 

It is known in the case of the Px +ipy superconductor in 
2D that a nontrivial topological charge of the free space 
Hamiltonian has a close connection to the existence of a 
localized zero energy state in the presence of a vortex [16, 
28]. We will see in the subsequent sections that the same 
correspondence is true in our Hamiltonian (3) describing 
superconducting Dirac fermions in 3D. We shall work in 
the chiral representation: 



and 



7°7 



cr 
-a- 



7 



/3=7° 



1 

1 



1 
-1 



(10) 



(11) 



III. EVEN PARITY PAIRING 

We first consider the case where the pairing takes place 
in the even parity channel, which is relevant to the color 
superconductivity of quarks [29, 30] . 



A. Topological chcirge of a free space Hamiltonian 

When the pairing gap is a constant, we can choose it 

to be real; A(a;) = A*(a;) = Aq. From Eq. (3), the free 
space Hamiltonian in the momentum space is given by 



Hp — 



a - p + /3m 
Ao 



Ao 

-cx - p — Pm + fj, 



Its energy eigenvalues have the usual form 



Ep = ±-^(Vm2+p2±^) +A2 



(12) 



(13) 



and each of them are doubly degenerate (signs arc not 
correlated). Note that the spectrum is fully gapped as 
long as Ao ^ 0. The computation of its topological 
charge is lengthy but straightforward.^ In order to eluci- 



date the effect of the fcrmion mass, we shall present the 
results for the m = case and the m ^ case separately. 



1. Chiral limit m = 

In the chiral limit m = 0, because the right-handed 
sector and the left-handed sector of the Hamiltonian (12) 
are decoupled, the unitary matrix qp in Eq. (8) has the 
block diagonal form: 



(14) 



Accordingly, we can define the topological charges for 
the right-handed sector and for the left-handed sector 
independently. The results are* 



and 



1 



247r2 
Ao 
2|Ao| 



dpe'^'^Ti 



247r2 

Ao 



dpe'^'^Tv 



2|Ao| 



(QRp^iqRp) i<lRpdjqRp)iqRpdkqRp) 

(15) 

(qlpdiqLp) {qZpdjqLp) {qlpdkqLp) 

(16) 



We find that each sector is topologically nontrivial having 
the nonzero topological charge.^ This implies the exis- 
tence of a localized zero energy state for each sector in 
the presence of a vortex. However, their signs are oppo- 
site and the total topological charge of the Hamiltonian 
(12) is vanishing; N = Nr + Nl = 0. 



2. Nonzero fermion mass m ^ 

When the fcrmion mass is nonzero m 7^ 0, the right- 
handed and left-handed sectors are coupled and thus the 
only total topological charge N is well defined. Because 
the spectrum is fully gapped for Ao 7^ 0, the inclusion of 
the fermion mass can not change the topological charge. 
Therefore we find 



N = 



(17) 



for arbitrary m, which means that the system is topolog- 
ically trivial. In the following subsection, we will see how 
these observations in the free space are reflected in the 
spectrum of fermions localized around a vortex line. 



* The topological charge for the full Hamiltonian (1) can be ob- 
^ Note that the topological charge is invariant as long as the spec- tained by simply summing contributions from all gapped sectors, 

trum is fully gapped. Therefore one can set n —> to reduce the ' The half-integer value of is common to relativistic fermions 
computational complication. because qp is noncompact at \p\ — >■ 00. See, e.g., Ref. [20]. 
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FIG. 1: Low-energy spectrum of fermions localized around the vortex line in the chiral limit (left panel) and with the nonzero 
fermion mass for the even parity pairing (right panel). The E < part of the spectrum is redundant in the superconductor 
and thus shown by the dashed line. 



B. Spectrum of fermions localized around a vortex 

line 



The spectrum of fermions in the presence of a 
singly quantized vortex line is obtained by solving the 
Bogoliubov-de Genne equation: 



—ia ■ d + j3m — fj, 
e-'''|A(r)| 



e*^|A(r)| 
ia ■ d — j3m + /i 



^{x) ^ E^{x). 

(18) 

Here we assumed that the vortex line extends in the z di- 
rection and A(a;) does not depend on z\ A(a;) = e*^|A(r)| 
where (r, 9, z) are cylindrical polar coordinates. Note 
that we do not make any assumptions on the form of 
|A(r)| except that it has a nonvanishing asymptotic 
value; |A(r — >■ oo)| > 0. Therefore the existence of lo- 
calized fermions that we will find below is independent 
of the vortex profile and thus in this sense they are uni- 
versal. This would be because these solutions have topo- 
logical origins and, in particular, the zero energy solu- 
tions are guaranteed by the index theorem [18, 31]. In 
contrast, there will be other Caroli-de Gennes-Matricon- 
type bound fermions on the vortex line which typically 
have the energy gap ^ |A(oo)p//i [32]. Because their 
spectrum depends on the vortex profile, we shall not in- 
vestigate such nonuniversal solutions in this paper. 

Because of the translational invariance in the z direc- 
tion, we look for solutions of the form 



<i>(r,0,z) = e'f^^0p>,0). (19) 
We rewrite the Hamiltonian in Eq. (18) as 



1. Zero energy solutions at m = ~ 

Consider the zero energy Bogoliubov-de Genne equa- 
tion a,t m = Pz = 0; Ho^o = 0. We can find two expo- 
nentially localized solutions (see Appendix A): 



ie^^Jiilir) 


e'7 

e-"'Ji(H 
~iJo{^j.r) 





-/J-|A(r')|dr' 



(21) 



and 



e 4 



g-/J-lA(r')|dr' 






-ie*Vi(^r) 



e-^^Ji(MO 



where A is a normalization constant: 

/•OO 

X = 2n drr [2J^{nr) + 2J^{fir)] e-^/o" |A(r-')l<ir' 
Jo 



(22) 



These two solutions have definite chirality; ^^(t>R/L 
i0i?/Li and hence their index. 



(23) 



ctzPz + fim 

-ctzPz - /3m 

(20) 



We first construct zero energy solutions for "Hq at m = 
Pz = and then determine their dispersion relations with 
treating SH {m,pz 7^ 0) as a perturbation. 



2. Perturbations in terms of m and pz 

We now evaluate matrix elements of SH with respect 
to 4>fi and It is easy to find 



27r foc 



JO 



Pz —im 
im —pz 
(24) 
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FIG. 2: The velocity v in Eq. (27) as a function of A//i. 



where we defined the parameter dependent velocity l^l < 
1 in units of the speed of hght: 



j^drr [4{^ir) - Jl{^lr)\ e-^/o |A(r')l<ir' 
j^drr[Jl{^Jir) + J2(^r)] e-^/o" \Hr')\dr- ' 



(25) 



Therefore, when to — 0, the right-handed and left- 
handed ferniions localized around the vortex line have 
the gapless dispersion relations: 



E ■ 



VPz 



and 



E 



(26) 



respectively. They have opposite velocities and thus 
propagate in opposite directions along the vortex line 
(Fig. 1). In the simple case where |A(r)| = A > is 
a constant, the velocity v in Eq. (25) can be evaluated as 



+ A2 £(-^2/^2) - /C(-/i2/A2) 



- 1 



A ^0 



(27) 



1 



3 

4 VA^ 



if 



oo 



which is plotted in Fig. 2 as a function of A//i. Here K, {£) 
is the complete elliptic integral of the first (second) kind. 

On the other hand, when to 7^ 0, the right-handed 
and left-handed fermions are mixed and their spectrum 
exhibits the mass gap provided by vm (Fig. 1): 



E 



Pi 



(28) 



This is closely linked to our previous observations on the 
topological charge of the free space Hamiltonian. In the 
chiral limit to = 0, each of the right-handed and left- 
handed sectors is topologically nontrivial {Nn — —Nl ^ 
0) and thus the vortex line supports the localized gapless 
fermions. However, once the fermion mass m 7^ is 
introduced, the total Hamiltonian is topologically trivial 
{N = 0) and thus the vortex line no longer supports the 
gapless fermions. Nevertheless, as long as vm <^ A(r — >■ 
c») is satisfied, the mass gap of such localized fermions 
is much smaller than the energy gap of bulk fermions 
and thus they can be important low-energy degrees of 
freedom. 



3. Nonperturbative solutions at fi = 

The above results rely on the perturbations in terms of 
TO and Pz- In the special case where = 0, we can obtain 
the exact dispersion relations of the localized fermions 
and their eigenfunctions with arbitrary to and Pz (see 
Appendix A). The two solutions to the Bogoliubov-de 
Genne equation (18) at fi = are found to be 

/ Pz+E± \ 




V 



TO 




-i{Pz + E±) 


«TO 



pipzZ-fg \A(r')\dr' 



(29) 



with E± = ± to2 + p2 _ ^jj^g discussions given above are 
valid in this special case too. 



C. Effective ID Hamiltonian along a vortex line 

Because the bulk fermions are gapped, the low-energy 
effective Hamiltonian of the system in the chiral limit 
m — should involve the gapless fermions existing along 
the vortex line. In order to write down the effective 
ID Hamiltonian, we expand the fermion operator in 
Eq. (3) in terms of the eigenfunctions of H: 

' dpz 
27r 



(30) 



where and bp^ are quasiparticle operators associated 
with the gapless right-handed and left-handed fermions, 
respectively. Because of the pseudoreality condition ^E" = 
C\l/* and the property of the solutions C0^y^ = (f>R/L, 

we have aj^^ = a_p^ and bj,^ — b^p^ . The quasiparticle 
operators obeying such conditions are called as Majo- 
rana fermions in condensed matter literatures [16, 19- 
22]. From the Hamiltonian H = ^ Jdx ^E'^'H^', the effec- 
tive ID Hamiltonian becomes 



V [dpz , I 
= 2 / 2^ (P^^P.^p. 



■Pzhlbp^ 



(31) 



When the nonzero fermion mass to 7^ is present, there 
are additional terms mixing a^^ and bp^ : 



Hid = \ {Pzal^ttp^ - Pzblbp^ 



+imb^p^ ttp^ — ima)p^ bp^ ) . 



(32) 



IV. ODD PARITY PAIRING 

So far we have discussed the case where the pairing 
takes place in the even parity channel. The situation is 
qualitatively different when the pairing takes place in the 
odd parity channel [A(a;) — ^ 7^A(a;)]. 
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A. Topological charge of a free space Hamiltonian 

When the pairing gap is a constant A(a;) = A* (a;) = 
Ap, the free space Hamiltonian in the momentum space 
is given by 



a ■ p + j3m — II 



T'^Ao 



7^Ao — a-p — /3m + /x 
Its energy eigenvalues have the form 



(33) 



± 



2 + p2 + ^2 + A2 ± 2^m2 (/x2 + A2) + p2^2 

(34) 

and each of them arc doubly degenerate (signs arc not 
correlated). Note that even if Aq ^ 0, the energy gap in 
the spectrum is closed at = ^i"^ + Aq. This fact will 
have important consequences on the topological charge 
of the free space Hamiltonian and thus on the existence 
of localized zero energy states in the presence of a vortex. 



1. Ghiral limit m = 

In the chiral limit m = 0, we can define the topo- 
logical charges for the right-handed sector and for the 
left-handed sector independently. Compared to the even 
parity pairing, 7''' in front of Aq flips the sign of the pair- 
ing gap only in the left-handed sector. Therefore, from 
Eqs. (15) and (16), we easily obtain 



Nr = 



Aq 
2|An 



and 



Nr, = 



Aq 

2|Ao| 



(35) 



We find that each sector is topologically nontrivial having 
the nonzero topological charge. The striking difference of 
the odd parity pairing from the even parity pairing is that 
their signs are the same and thus the total topological 
charge of the Hamiltonian (33) is nonvanishing; 



N = Nr + Nl = 



An 



^0 



(36) 



This implies, unlike the even parity pairing, that the ex- 
istence of localized zero energy states in the presence of a 
vortex is robust against the inclusion of the small fermion 
mass. 



another fully gapped phase characterized by its vanishing 
topological charge^: 



N : 



for w? < tJ? + A^ 



- (37) 
for m2 > jji^ + Aq. 



Because these two phases can not be distinguished by 
symmetries, the phase transition between them is a topo- 
logical phase transition. This is exactly the same type 
of the topological phase transition occurring in the 2D 
Px + iPy superconductor at ^ = as a function of the 
chemical potential [16, 28]. Other quantum phase tran- 
sitions resulting from the momentum space topology are 
extensively discussed in Rcf. [33]. In the following sub- 
section, we will see how these observations in the free 
space are reflected in the existence of gapless fermions 
localized around a vortex line. 



B. Spectrum of fermions localized around a vortex 

line 



The spectrum of fermions in the presence of a 
singly quantized vortex line is obtained by solving the 
Bogoliubov-de Genne equation: 

•fe "^|A(r)| lOL - d - l3m-\- jij ^ ' ^ ' 

(38) 

As we stated in Sec. IIIB, we do not make any assump- 
tions on the form of |A(r)| except that it has a nonva- 
nishing asymptotic value; |A(r — >■ 00 )| > 0. Because of 
the translational invariance in the z direction, we look 
for solutions of the form 

^{r,e,z) = e'P^'(i>j,Xr,e). (39) 

We rewrite the Hamiltonian in Eq. (38) as 



-'P'^ne'p^' = n\p^^Q 



OtzPz 

-a^p, 



(40) 



= Ho + 5H. 



We first construct zero energy solutions for Hq at = 
and then determine their dispersion relations with treat- 
ing 5% {pz 7^ 0) as a perturbation. 



2. Nonzero fermion mass m =^ 



Zero energy solutions at pz = 



When the fermion mass is nonzero m ^ 0, the only to- 
tal topological charge N is well defined. Because the 

spectrum remains gapped as long as rri^ < fi^ + Aq, 
the topological charge remains the same; N = Ao/|Ao|. 
However, when the fermion mass exceeds the critical 
value to2 = ^2 ^ ^2 g^-^ which the energy gap in the 
spectrum closes, there is a quantum phase transition to 



Consider the zero energy Bogoliubov-de Genne equa- 
tion a.t Pz = 0; Ho<po = 0- We can find two potentially 



This phase is continuously connected to the Hamiltonian with 
/I, Ao — >^ where the topological charge can be computed most 
easily. 
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normalizable solutions (see Appendix B): 



2. Perturbation in terms of pz 



e 4 

7x 



and 



e 4 

7\ 



is/ IJp' — m? e*^ Ji 

0_ 

V'm^ - m2_e-*^ Ji 

_ 

\ -imJo / 



mJo 
0_ 



g-/J-|A(r')|dr' (4;^) 



V 



_ 

— imJo 



g-/J-|A(.')l<ir'^ (42) 



where we introduced shorthand notations; Jq = 
Jo(-\/ /i^ — r) and Ji = Ji(-\/ yu^ — r), and A is a 
normalization constant: 



/■OO 

A = 27r / (ir r [2(/i2 + m^) + 2(;U^ - m^) J?] 

^0 



X e 



-2/J-|A(r')|c(r' 



(43) 



These two solutions in the chiral limit m — i) have definite 
chirality; "f^(t>R/L — i'^fi/i) ^nd hence their index. Note 
that, although the nonzero fermion mass m ^ mixes the 
right-handed and left-handed fermions, their gaplessness 
is preserved. 

We now examine if the above two solutions are normal- 
izable or not. When rr? < fi^, they are normalizable ow- 
ing to the exponentially decaying factor ■I'o 1^^'' ^■ 



g-r|A(oo)| 

Jo = -fo(v^^ ~ /U^r) and Ji 



00. When m > ji , we note that 
il-i_{y/m? — fj,^ r) expo- 



nentially grow as Jo,Ji ^^/rn^-t^^r^ Nevertheless, as 
long as < fj,"^ + |A(oo)p is satisfied, the exponen- 
tially decaying factor dominates and the solutions are 
still normalizable. However, when the fermion mass ex- 
ceeds the critical value = /U^ -t- |A(oo)p, the two so- 
lutions exponentially grow and thus they are no longer 
acceptable solutions. Therefore we find that the gap- 
less fermions localized around the vortex line exist only 
when < n"^ + |A(oo)p and they disappear when 
> /z^ -)- |A(oo)p. 
Of course this is closely linked to our previous ob- 
servations on tlic topological charge of the free space 
Hamiltonian; it is topologically nontrivial (iV ^ 0) for 



|Aor and trivial (iV = 0) for w/' > + 



rri^ < p'^ - 

and there is a topological phase transition in between. 
Also note the striking difference of the odd parity pair- 
ing from the even parity pairing where the presence of 
the fermion mass immediately opened up a mass gap for 
the fermions localized around the vortex line. 



We now evaluate matrix elements of 5% with respect 
to and (j)L for m? < fj? + |A(oo)p. It is easy to find 



27r 

de 

JO 



Pz 

-p. 



V = 



(44) 

where we defined the parameter dependent velocity \v\ < 
1 in units of the speed of light: 

J^drr [ifi^ + m^)J§ + {fi^ - m^)Ji] e-^SS \Mr')\dr' ' 

(45) 

Therefore the dispersion relations of the gapless fermions 
localized around the vortex line are given by 

E = vpz and E = —vpz- (46) 

They have opposite velocities and thus propagate in op- 
posite directions along the vortex line. In the simple 
case where |A(r)| = A > is a constant, the velocity v 
in Eq. (45) can be evaluated as 

V = ^ ^ ^ - 1 (47) 

(/i2 + A2)5-(/z2 + A2-m2)/C ' ^ ' 

where K = IC{{m^ - ^^^)/A^) and £ = £{{m^ - ^i^)/A^). 
When m — 0, this is identical to Eq. (27) and plotted 
in Fig. 2 as a function of A/fi. The velocity vanishes at 
m = fi and changes its sign for m > p. 



3. Nonperturbative solutions at n = 

The above results rely on the perturbation in terms of 
Pz- In the special case where /x = 0, we can obtain the 
exact dispersion relations of the localized fermions and 
their eigenfunctions with arbitrary pz (see Appendix B). 
The two solutions to the Bogoliubov-de Genne equation 
(38) at = are found to be 

/ Io{mr) \ 



e'T ze'^7i(mr) 

7^ 

— i/o(mr) 
e~*^ Ji (mr) 

V / 



giPzZ-J^ \A{r')\dr' j-^g-j 



with E = Pz and 



^{r,0,z) 



7^ 



( \ 

-ie Ii{mr) 
Io{mr) 


-e~^^Ii{mr] 



-iIo{mr) ) 



JPzZ-Jg \A{r')\dr' 



(49) 
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with E = —pz- They are normahzable in the x-y plane 
as long as m < A(oo). The discussions given above are 
valid in this special case too. 



C. Effective ID Hamiltonian along a vortex line 



Because the bulk fermions are gapped, the low-energy 
effective Hamiltonian of the system with < /x^ + 
|A(oo)p should involve the gaplcss fermions existing 
along the vortex line. In order to write down the ef- 
fective ID Hamiltonian, we expand the fermion operator 
\E' in Eq. (3) in terms of the eigenfunctions of H: 



where a„ and 6„ 



(50) 



are Majorana quasiparticlc operators 
. and h]p__ = b-p^ and associated with 



obeying a],^ = a_p 

the two gapless fermions. From the Hamiltonian H 
^ Jdx 'i/^'H'^, the effective ID Hamiltonian becomes 



BOUNDARY PROBLEMS 



(51) 



1. Chiral limit m = 

In the chiral limit m = 0, we can find gapless right- 
handed and left-handed fermions localized at the bound- 
ary. Their dispersion relations are given by 



E = ± 



(53) 



where p±_ = {px,Py) is the momentum parallel to the 

boundary and |?,'| < 1 is the parameter dependent veloc- 
ity in units of the speed of light: 



V 



I ^2i^LZ-2 Aiz')dz' 

J —oo 

J — OO 



(54) 



In the simple case with A{z) = Aoo sgn(z), the velocity 
can be evaluated as u = A^/ (A^ -h/i^). The low-energy 
effective Hamiltonian involving the two gapless fermions 
existing at the 2D boundary becomes 



H2D = Y 



bl, (o-_L ■p±)bp 



with <T 



t 



Here a„ , and bn 



(55) 
arc two- 



, (J2). j._n^±i^ "■p± "ii"-! "p± 

component Majorana quasiparticlc operators obeying 

6T 



a. 



p± 



al_„ . cTi and 6j, 



-p± 



ax. 



In Sees. Ill and IV, we have discussed the connec- 
tion between the nonzero topological charge of the free 
space Hamiltonian and the existence of localized gapless 
fermions in the presence of a vortex. Another character- 
istic of the topological state of matter is the existence of 
gapless edge/surface states when it is terminated by an- 
other gapped state having a different topological charge. 
In this section, we will investigate two types of boundary 
problems and elucidate the existence of gapless surface 
fermions localized at the boundary. 



A. Boundary between Ao > and Ao < 

We first consider the boundary where the pairing gap 
A(x) = A(^) changes its sign as (A^o > 0): 



A(^)^ 



-l-Aoo for ^ 00 
—Aoo for z —00. 



(52) 



Because the topological charge of the free space Hamilto- 
nian, when it is nonzero, depends on the sign of the pair- 
ing gap [see Eqs. (15), (16), and (35)], the above bound- 
ary connects two gapped phases with different topological 
charges and thus we expect the gapless fermions localized 
at the boundary. Here we present the low-energy spec- 
trum of such fermions and their effective 2D Hamiltonian. 
Their derivations arc analogous to those in the vortex 
problems and some details are provided in Appendix C. 



2. Even parity pairing with m ^ 

As we can expect from the topological charge of the 
free space Hamiltonian [Eqs. (17) and (37)], the effect of 
the nonzero fermion mass to 7^ is qualitatively different 
for the even parity pairing and the odd parity pairing. 
When the pairing gap is parity even, the presence of the 
fermion mass immediately opens up a mass gap for the 
fermions localized at the boundary and the dispersion 
relations become 



E = ±\v\ 



+pI+pI- 



(56) 



This is linked to the fact that the two phases at both 
sides z ±00 have the same topological charge N = 0. 
Now the effective 2D Hamiltonian has additional terms 
mixing Op^^ and br. 



^p± ■ 



H. 



2D 



3. Odd parity pairing with m ^ 



(57) 



On the other hand, when the pairing gap is parity odd, 

the two phases at both sides have different topological 
charges TV = ± 1 even in the presence of the small fermion 
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mass. Accordingly, the boundary still supports the two 
localized gapless fermions [Eqs. (53) and (55)] with the 
modified velocity < 1: 



for 



< /Lt^ and 



J — oo 



V = 



m 



(58) 



(59) 



for /i^ < < /Lt^ + A^. 

When the fermion mass is increased further and ex- 
ceeds the critical value = /i^ + A^, the resulting 
solutions are no longer normalizable [sec Eqs. (C3) and 
(C9)] and the gapless surface fermions disappear. This 
is linked to the fact that the two phases at both sides 
z ±oo have the same topological charge A'' = as a 
consequence of the topological phase transition. 



B. Boundary between superconductor and vacuum 

We then consider the following boundary that models 

the interface between the superconductor at z < and 
the chiral symmetry broken vacuum at z > 0: 



/u, A > 0, TO = for 2 < 
^ = A = 0, TO > for z > 0. 



(60) 



For the even parity pairing, the two gapped phases at 
both sides have the same topological charge N = 0, while 
they have different topological charges for the odd parity 
pairing; A^ = 1 at z < and A^ = at z > 0. 

Accordingly, when the pairing gap is parity odd, we 
can find one gapless fermion localized at the boundary 
z = 0. The corresponding two zero energy solutions at 
Px=Py = are given by 







i ^ \ 












































V ; 







(61) 



for z < and 



4>(a;) oc 














1 


—i 












i 


i 












—i 
















(62) 



for z > 0. Its effective 2D Hamiltonian becomes 



2D 



\V\ f dpj_ r . -, 



(63) 



1 



/(I 



and the Ma- 



' A+ifj, 

aE.p^ai. However, when the 



with the velocity v = 
i or ana condition al, 

pairing gap is parity even, we cannot find such localized 
zero energy solutions that are continuous at the bound- 
ary z = and therefore gapless surface fermions do not 
exist. 



VI. SUMMARY AND CONCLUDING 
REMARKS 

Motivated by the recent discoveries of the topological 
insulators and superconductors, we have investigated the 
topological properties of superconducting Dirac fermions 
in 3D both for the even parity pairing and the odd parity 
pairing. The results are summarized in Table I. In the 
chiral limit to = 0, we find that the system is topologi- 
cally nontrivial in the sense that each of the right-handed 
and left-handed sectors of the free space Hamiltonian has 
the nonzero topological charge; Nji^L 7^ 0. Accordingly, a 
vortex line in the superconductor supports localized gap- 
less right-handed and left-handed fermions. Their dis- 
persion relations are given hy E = ±vpz, where |u| < 1 
defined in Eq. (25) is the parameter dependent velocity, 
and thus they propagate in opposite directions along the 
vortex line. 

The effect of the nonzero fermion mass to ^ is quali- 
tatively different for the even parity pairing and the odd 

parity pairing. When the pairing gap is parity even, the 
presence of the fermion mass immediately opens up a 
mass gap for the fermions localized around the v ortex line 
and the dispersion relations become E = ±vy^m? + p^. 
This can be understood from the vanishing total topo- 
logical charge of the free space Hamiltonian for the even 
parity pairing; N{= Nr + Nl) = 0. 

On the other hand, when the pairing gap is parity odd, 
the total topological charge of the free space Hamilto- 
nian is nonvanishing N and thus the system remains 
topological even in the presence of the small fermion 
mass. Accordingly, the vortex line still supports the lo- 
calized gapless fermions. When the fermion mass is in- 
creased further, there is a topological phase transition at 
TO = + A^, where the topological charge jumps from 
the nonzero value to zero and consequently the gapless 
fermions localized around the vortex line disappear. 

Our results for the even parity pairing are relevant 
to the color superconductivity of quarks. In the CFL 
phase where all nine quarks are gapped, the mean-field 
model Hamiltonian for the color superconductor (1) has 
the topological charge Nr = —Nl = Y^j=i 2\i\-\ 
chiral limit to = 0. Therefore the U(1)b vortex line that 
arises when the CFL quark matter is rotated [34, 35] 
supports nine sets of localized gapless right-handed and 
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TABLE I: Summary of the topological charge of the free space Hamiltonian and the low-energy spectrum of fermions locahzed 
around a vortex line (or a boundary where the pairing gap changes its sign), u is a parameter dependent velocity. 





even parity pairing 


odd parity pairing 




topological charge mid-gap state dispersion 


topological charg 


e mid-gap state dispersion 


m = 






1 E = ±vpz 




A^ = E = ±v^m'^ +p'i 


\ AT = 


E = ±vpz (m^ < + Ao) 
none (m^ > fi^ + Aq) 



left-handed quarks. In the presence of the small quark 
mass or chiral condensate, such localized quarks become 
gapped but their mass gap vm is much smaller than 
the energy gap of bulk quarks A.'^ Furthermore, the 
mass gap at high density vm m (A/;u)^ ln(/z/A) [see 
Eq. (27)] is parametrically smaller than the masses of 
pscudo-Nambu-Goldstonc bosons ^ m (A//i) in the CFL 
phase [36], which are important to the transport prop- 
erties and neutrino emissivity of the CFL quark mat- 
ter [24]. Whether such new low-energy degrees of free- 
dom localized around the vortex line have some impact 
on the physics of rotating neutron/quark stars is an im- 
portant problem and should be investigated in a future 
work. 

Our results for the odd parity pairing might be irrele- 
vant to the color superconductivity of quarks. Neverthe- 
less it would be interesting to consider if the topological 
phase transition found in this paper can be realized in 
condensed matter systems where the 2D Dirac fermions 
appear. 

Finally, Table I reveals the intriguing connection 
between the nonzero topological charge of the free 
space Hamiltonian and the existence of localized gap- 
less fermions in the presence of a vortex. We also eluci- 
dated the existence of gapless surface fermions localized 
at a boundary when two phases with different topological 
charges are connected. The mathematical proof of these 
correspondences remains an open question. 

Note added. — After the submission of this paper, 
there appeared a paper by S. Yasui, K. Itakura, and 
M. Nitta [37] which has some overlap with the present 
paper. In their paper, a non-Abelian vortex in the CFL 
phase is also considered. 
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Supposing ~ 500 MeV, A ~ 50 MeV, and m ~ 10 MeV and 
using Eq. (27), we can estimate the mass gap to be as small as 
vm/A ~ 5 X 10~^. 



Appendix A: Derivations of solutions for even parity 
pairing 

Here we outline how the solutions to the Bogoliubov- 
de Genne equation (18) for the even parity pairing are 
derived. We introduce notations 



Fl 
Gr 
\GlJ 



(Al) 



where Fjnj^) and Gm^^ arc right-handed (left-handed) 
two-component fields. We first make an ansatz 

Fn{r, 6) = [js^'l^^r)) ^"^'^ 



and 



Gfl(r, e)=r^ g- /; |A(.')|rfr' (A3) 

V 9Ri{r) J 

and the same fov R ^ L so that they are exponen- 
tially localized in the x-y plane. Then we look for /^(l) 
and gii(L) that are regular at origin and independent of 
]A(r)|. |A(r)| can be eliminated from the equations by 
imposing 



9r 



-ic^ifn and gL=i<TifL- 



(A4) 



Now fa and /l satisfy the following four sets of equa- 
tions: 



(A5a) 
(A5b) 
(A5c) 
(A5d) 








E 



fL + mfR = 0. 



We can find consistent solutions in two cases; m = Pz = 
E = [Eqs. (21) and (22)] and /x = [Eq. (29)]. 
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Appendix B: Derivations of solutions for odd parity 
pairing 

Here wc outline how the solutions to the Bogoliubov- 
de Genne equation (38) for the odd parity pairing are 
derived. We use the same notations as in Appendix A. 
In this case, |A(r)| can be eliminated from the equations 
by imposing 

gn = -iaifR and gL = -icrifb- (Bl) 

Now fn and /i, satisfy the following four sets of equa- 
tions: 



Pz-E 
, -Pz-E, 



-Pz-E 
Pz-E, 



i = (B2a) 

/ij = (B2b) 

= (B2c) 

/l = 0. (B2d) 



We can find consistent solutions in two cases; Pz = E = 
[Eqs. (41) and (42)] and /i = [Eqs. (48) and (49)]. 



Appendix C: Derivations of solutions for boundary- 
problems 

Here we outline how the spectrum of fermions in the 
presence of the boundary (52) studied in Sec. V A is ob- 
tained. The Bogoliubov-de Genne equation to be solved 
is 



-ia-d + pm — ^ (75)A(z) 



(75)A(z) 



iot-d — j3m + jjL 



$(x) = E^{x). 

(CI) 

Because of the translational invariance in the x-y plane, 
we look for solutions of the form 



$(x) = e^^''-+'^yy<t>j>^,p^{z). 
We first make an ansatz 



(C2) 



Il{z) 
9r{z) 

\gL{z)J 



- A{z')dz' 



(C3) 



so that the solution is exponentially localized in the z 
direction. Then we look for fn^L) and gii(L) that are 
independent of /S.{z). /^{z) can be eliminated from the 
equations by imposing 



where the upper (lower) sign corresponds to the even 
(odd) parity pairing. Now and satisfy the following 
four sets of equations: 

fR + m/i = (C5a) 
fn ±mfL=0 (C5b) 
/l + m/fl = (C5c) 
fL ± m/K = (C5d) 

T + M - 
with p± =px±ipy. 

1. Even parity pairing 

When the pairing gap is parity even (upper sign), we 
can find the following four zero energy solutions at m = 
Px=Py = 0: 





E 




p- 


p+ 


— 


i 


- - 




-E 


i 


P- 


p+ 




+ M- 






i 


~P- 


-p+ 




-fj,- 




E 




-P- 


-p+ 




§^ 

i 


+ M- 



oc 







( ' \ 


( " \ 


( ' \ 

























' 


V ) 




I J 


V ; 





(C6) 



By evaluating matrix elements of the perturbation 
Hamiltonian 

SH= ( "^^^ "^^^ + /3m \ 

y -a^Px - Q-yPy - pmj 

(C7) 

with respect to the corresponding four zero energy eigen- 
functions (^0,0(2) in Eq. (C3), we obtain the following 
effective 2D Hamiltonian: 



H2D 



/ v*p- 
vp^ 
vm 

\ v*m 



v*m \ 

vm 

—vp. 

-v*p+ / 



(C8) 



gR 



-icr^fR 



and 



gL = iiaa/i, 



(C4) 



with the velocity v defined in Eq. (54). Its energy eigen- 
values are given hy E = ±.\v\^rTn? +p%+Py. This spec- 
trum with ti = 1 becomes exact for arbitrary Px,Py, and 
m in the special case where = 0. 



2. Odd parity pairing 

When the pairing gap is parity odd (lower sign), we 
can find the following four zero energy solutions at Px = 



12 



Py=0: 







V 



m 






IJ.+ - 


V m / 



m 




\ 
I 






m 




\lJb + ^/Jjfl^^rr? ) 



\ 



J 



(C9) 



Note that the corresponding four zero energy cigenfunc- 
tions (/)o.o(2) in Eq. (C3) are normalizablc in the z direc- 
tion as long as < /x^ + A^. By evaluating matrix 
elements of the perturbation Hamiltonian 



5H = 



axPx 





~^xPx ~ ^yPy ) 



(CIO) 



with respect to the four zero energy eigenfunctions, we 
obtain the following effective 2D Hamiltonian: 



/ v*v. 

vp+ 



V 








-v*p+ 










(Cll) 



for m < /X with the velocity v defined in Eq. (58) and 



U2B 



( 








—ivp_ 


A 








—ivp+ 










ivp- 










\ivp+ 














(C12) 



for n"^ < < /U^ + with the velocity v defined in 
Eq. (59). Their energy eigenvalues are given by £^ = 

±.\v\^p% +Py- This spectrum with v = 1 becomes exact 

for arbitrary p^ and Py in the special case where /U = 0. 
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